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ABSTRACT 


\ 
A new approach to defining inlet distortion factors for axial-flow 


compressors is presented. A summary of past analyses of inlet distortion, 
along with a brief review of flow-distortion definition and compressor 
performance, provides a background for the proposed theory. Inlet flow- 
distortion in terms of total-pressure fluctuations at the compressor face 
is converted to vorticity. The effects of vorticity on the internal 
aerodynamics of the compressor are investigated with some approximate 
calculations included to indicate the validity of this approach. Refine- 
ments and alternate solutions to the theory which could lead to definition 
of a more comprehensive and reliable definition of stall-inducing inlet 


distortion are included. 
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I. INTRODUCTION 


Engine stall has always been a problem of varying magnitude in jet 
powered aircraft employing axial-flow compressors. At present aircraft 
gas turbine performance is largely limited by flow separation and surge 
in the compressor. Generally, highest compressor efficiencies occur 
near the point of flow breakdown (stall), and so the design margin 
between normal operation and stall is kept to a practical minimum. A 
compressor's ability to remain below this stall limit is significantly 
affected by the amount and degree of distortion in its inlet flow. 

Inlet flow distortion has been called by many names, all of which basi- 
cally describe variations in total-pressure, total-temperature, or 
velocity over the inlet cross-section. Previously, inlet flow distor- 
tion has been at a level which the compressor could handle without a 
serious degradation in performance. However, with the greater speeds 

and performance demanded of today's high-specific-energy turbomachines, 
inlet distortion at the compressor face can alone be great enough to 

cause compressor stall and surge. One prominent source of distortion 

has been traced to shock-wave/boundary-layer interactions within the 
inlet. Others are steam ingestion, pressure transients due to afterburner 
light-off, ordnance discharges, etc. The magnitude of the problem in 

some recent aircraft designs has caused much research to be conducted in 
this area, specifically correlating maps of total-pressure variations at 
the compressor face to engine surge. The results of this work have 
yielded an immense amount of empirical data. However, reduction of these 
data has yet to provide a technique which adequately describes the distor- 


tion accommodation of a given engine-inlet configuration over its entire 
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range of operation. In contrast to empirical correlations yielding a 
Statistical yes-no answer to the question of whether or not the engine 
will stall, methods are needed which link stall criteria for a blade 
element to distortion patterns. This paper presents a new approach to 
analyzing already existing data and provides evidence of a better 


insight into the actual fluid mechanics involved in the problem. 
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II. STATEMENT OF PROBLEM 


A. AXIAL COMPRESSOR PERFORMANCE! 

The performance of an axial-flow compressor such as found on air- 
craft turbojet engines is best described by its steady-state performance 
map. Figure 1 illustrates a typical compressor map. "Pressure ratio" 
is exit total-pressure to inlet total-pressure and "corrected airflow" 
is engine mass-flow rate corrected for total-temperature and total- 


pressure. 


Surge Line 







Characteristics 









Nominal 


AAA —- 


Operation 


Constant Speed 
Lines 


Stall /Surge Margin 


COMPRESSOR PRESSURE RATIO 


l4 S 16 17 8 19 20 2 
CORRECTED AIRFLOW (kg/sec) 


FIG.1 COMPRESSOR PERFORMANCE MAP 
J85- GE- I3 [adapted from Ref.2] 


IThe information and examples in this section were adapted from a 
report by Brimelow [Ref. 1]. 
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The surge line defines the maximum pressure ratios which the compressor 
can achieve throughout its range of rotor speeds. Operation above this 
line is accompanied by a drastic reduction in pressure ratio. This is 
called "surge," and it is caused by excessive aerodynamic blade stall. 
The stability of a system is determined by the slope of its pressure 


ratio/flow characteristic. Consider steady operation at point М, of 







are 2. 
Compressor Characteristic 
О 
E 
a 
u) 
аа 
Ф Receiver p 
ul Characteristic 
а 
CORRECTED AIRFLOW 
FIG. 2 COMPRESSOR / RECEIVER 


MATCHING CHARACTERISTICS 


Now if the flow became unsteady and a perturbation were to cause a 
reduction in airflow to point S4 (constant speed), the compressor would 


be delivering an instantaneous pressure greater than that in the 
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receiver. Effectively the compressor back-pressure has been lowered, 

and the flow would increase until a match is again reached. Similarly 
for an increase in flow to point S5, the increased back-pressure would 
cause a return to stable operation at point M. However, if the compres- 
sor is matched on the positive slope of the characteristic, by the same 
reasoning, unsteady conditions exist and perturbations are aggravated 

to "surge." 

A stall margin, defined as the vertical distance between the nominal 
(design) operating line and the surge line, can be thought of as a 
measure of the compressor's ability to accommodate flow transients which 
instantaneously load the machine. This margin is dictated by overall 
engine-inlet design and is affected by many different phenomena. Unfor- 
tunately the point of maximum efficiency in a compressor is just prior 
to stall or just below the surge line on the performance map, and there- 
fore a trade-off between efficiency and surge margin is required. The 
expected flow transients which the compressor must be able to accommodate 
will affect the design efficiency. 

Flow transients initiated within the engine system are normally 
caused by power level (throttle) movement. Transients can also be 
initiated by flow distortion at the compressor face. Of the various 
flow distortion parameters , Brimelow [Ref. 1] states that inlet total- 
pressure and total-temperature have been identified as having the greatest 
effect on engine stability, which is to a large extent the compressor's 
ability to remain below the surge line. Total-pressure variations can 
be steady-state spatial distortion, time-varying spatial distortion, 
in-phase pulsations which are spatially uniform over the compressor face, 


or a combination of these. Low frequency spatially uniform (planar) 
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pressure fluctuations (on the order of 20 CPS) do not affect the com- 
pressor but do have a detrimental effect on the inlet and engine control 
systems. High frequency planar fluctuations (more than one per revolu- 
tion) affect compressor stability much like time-varying spatial 
distortion. 

Spatial distortions can be further segregated into radial and cir- 
cumferential inlet flow distortion. This separation ínto circumferential 
and radial distortion was arbitrary, somewhat motivated by geometric 
simplicity. A radial total-pressure distortion such as shown in Fig. 
3(a) produces very little variation in inlet static pressure. The net 
result is a low velocity in the low pressure region and a high velocity 


in the high pressure region. 





LOW PRESSURE 


HUB TIP 
(а) (b) 


FIG. 3 RADIAL PRESSURE DISTORTION 
[adapted from Ref. ] 
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In the low velocity region the blades will operate at higher incidences 
and higher total-pressure ratios. The net result will be a reduction in 
the compressor stall margin. 

For the circumferential type distortion shown in Fig. 4(a), the 
associated surge line reduction is a function of the basic stage-match 
of the compressor, the axial spacing between blade and vane rows, and 
the unsteady response of the -—. blades. Since the blade rows restrict 
circumferential redistribution of flow, each circumferential segment 


performs essentially independently. 


Undistorted 
Inlet Flow 





Low Pressure 





(a) |809 360° 
(b) 


FIG. 4 CIRCUMFERENTIAL DISTORTION 
[ adapted from Ref.!] 


All segments are required to deliver the same exit static pressure so 
that in the area of low inlet pressure there is a high pressure ratio 


and vice-versa. Correspondingly, the high pressure ratio forces a low 
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axial velocity in this area and a high velocity in the segment of high 
inlet pressure. These variations in velocity, produced less than one 
duct-diameter upstream of the compressor, result in static pressure 
fluctuations at the compressor face. This mechanism of varying veloci- 
ties (incidences) can induce surge when the average pressure ratio of 
the distorted flow is less than the surge pressure ratio for undistorted 
inlet flow. The effects of circumferential pressure pulses on a TF-30 
engine are presented in Ref. 3. It was found here that the amplitude of 
a pulse necessary to stall the engine varied inversely with pulse dura- 
tion and that the engine was most sensitive to distorted sectors of 180 
to 240 degrees. 

The distortions which are encountered in actual compressor applica- 
tions such as jet aircraft inlets are neither radial nor circumferential 
but a combination of the two. Furthermore the magnitude, location, 
shape, and duration of these pressure variations change continually. 

Je is virtually impossible to evaluate precisely the effect that these 
distortions have on engine stability because there are so many distor- 
tion patterns significant to each compressor surge. Therefore it has 
been necessary through experimental programs to develop generalized 
correlations between pattern characteristics and corresponding effects 

on compressor stability. Correlating the shape, location, spatial 
extent, and magnitude of the inlet distortion requires analyzing many 
radial and many circumferential patterns and relating the overall picture 
to engine stability. From this is generated a distortion factor which is 
a measure of the associated loss in stability-margin for that given 
propulsion system. Because compressor geometry and blading are not 
considered in arriving at this parameter, the result cannot be applied 


| 
to other systems. | 
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For a more detailed and complete analysis of this subject, consult 


Brimelow's excellent work in Ref. 1. 


B. PRESENT STATE OF INLET DISTORTION ANALYSES 

Until a few years ago, engine-inlet compatibility limits were success- 
fully determined by evaluating an engine's tolerance to steady-state 
distortion. Normally this was done with screens or similar obstructions 
placed in the inlet and the compressor face instrumented with low-response 
pressure probes. The engine was operated under all conditions of interest, 
and by comparing pressure data at the compressor face with the engine's 
performance, distortion factors were derived. It was known that unsteady 
flow was being produced by these test devices, but no attempt was made to 
evaluate its effect. However, when this technique was applied to the 
TF-30 turbofan engine during the F-111A flight test program, the data did 
not correlate well nor was the method found valid for other fan type 
engines at that time. An improved steady-state distortion factor, Kap: 
was developed, but it too fell short of providing the needed total evalu- 
ation of inlet distortion. It was becoming more evident that there was 
another significant variable involved, namely turbulence. Turbulence can 
be described as flow with random pressure fluctuations and can be quanti- 
fied as a fluctuation in RMS total-pressure, (АРрм)/ (РТ). With 
increased aircraft speeds, turbulence has taken on greater importance 
since there is a significant amount of turbulent flow associated with 
the boundary-layer/shock-wave interaction within supersonic inlets. It 
was felt that steady-state distortion indexes were still of value, but 
a correlation with turbulence was necessary. One such empirical relation 
showing a decrease in distortion tolerance, Ко» with increasing turbulent 


energy is depicted in Fig. 5. 
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| adapted from Ref. 4 | 


Turbulence has two primary effects on axial compressor performance. 
First, turbulent inlet flow causes a reduction in compressor surge line, 
its sensitivity being a function of turbulent frequency and amplitude. 
Second, a reduction in flow capacity accompanies inlet turbulence. These 
two phenomena were observed in tests performed by Brimelow and Plourde 
[Ref. 5] on a TF-30 turbofan engine. This all led to the development of 
new distortion indexes such as the Kp, ‚and Kom parameters of Ref. 6. 
However, lack of correlation caused Burcham and Hughes to raise doubt as 
to the usefulness of the conventional turbulence definition, (АРрм)/ (РТ). 
Bellman and Hughes [Ref. 7] found no consistent relationship at all 
between turbulence factors and surge, and in some cases surge was experi- 
enced as turbulence level was decreasing. Therefore it was concluded here 
that peak-to-peak pressure fluctuations at key areas of the compressor 


face are more important in surge prediction than are turbulence levels. 
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Also extensive experimental data have been obtained from a J85-GE-13 
turbojet placed in a supersonic wind-tunnel running at Mach numbers 
close to 2.5. By varying the characteristics of the engine inlet duct, 
a desired level of time variant flow distortion was produced through 
internal shock interactions [Ref. 2]. Calogeras, Burstadt and Coltrin 
found for this J85 configuration that time-averaged distortion could be 
well below a stall- inducing value while instantaneous peaks exceeded 
this value. In addition, severe pressure fluctuations could result in 
an instantaneous distortion of a type quite different from the time- 
averaged distortion. 

Of these few representative programs discussed, there are two 
striking similarities worthy of note: (1) Each is dedicated to using 
inlet flow distortion described by raw total-pressure data to predict 
compressor surge. No attempt at determining the cause or nature of the 
stall is made. (2) There is much uncertainty and disagreement associated 
with the present methods of predicting distortion induced stall. A 
universal distortion parameter has not been agreed upon. 

Obviously a better understanding of the internal fluid dynamics, 
while a difficult task at best, would nonetheless allow a more compre- 
hensive generation of a stall-inducing distortion parameter.  Total- 
pressure is the easiest and perhaps the only flow parameter which can be 
measured at the compressor face with sufficient accuracy and frequency 
response to gain meaningful data. The vast amount of pressure data 
available due to the previously mentioned programs provides an exten- 
Sive empirical data base with which to treat this problem. In 
evaluating these data, inclusion of time and temperature dependence 
could only improve the usefulness of the results. Therefore it is 


proposed that vorticity become the basis of a "universal distortion 
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factor." This thesis will show how existing total-pressure maps along 
with total-temperature and time data can be transformed into vorticity 
maps. Analysis of the resulting vorticity maps will provide evidence 
as to the flexibility of this method. 

Vorticity is more closely related to fluid dynamics than is total 
pressure alone. Compressor stall and engine surge are flow-separation 
phenomena; therefore the key to understanding their causes lies in the 


internal dynamical situation. 


C. VORTICITY APPROACH 
1. Compressor Stall 

The most significant phenomenon affecting the performance of an 
axial-flow compressor is stall. If the flow rate through a compressor 
blade row is decreased while the rotor speed is maintained constant, 
angles of attack on the blades will increase and eventually flow sepa- 
ration will occur. This condition is called stall and is similar in 
some respects to the stall of an isolated airfoil; differences are due 
to the adjacent blading in a cascade. Stalled operation of a compressor 
is accompanied by a drastic decrease in efficiency and excessive oscil- 
lating blade loads capable of causing structural failure. Stall occurs 
in the range of unstable compressor operation and hence can trigger 
Surge. Usually a compressor stalls in an asymmetric manner; that is, 
well defined regions of the annulus are stalled while the remaining area 
is unstalled; see Fig. 6. The flow through these stalled patches is 
severely restricted in accordance with the decrease in total flow-rate 
through the compressor, thus the flow is diverted around this blockage. 
Effectively this increases blade angle-of-attack on one side of the 


restriction while decreasing it on the other side. In Fig. 7 blade No. 1 
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FIG. 6 TYPICAL STALL PATTERNS 
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FIG. 7 STALL CELL GROWTH AND 
ROTATION IN A BLADE ROW 
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-will soon become unstalled while at blade No. 5 flow separation will 
commence. This condition is termed "rotating stall" since it propagates 
circumferentially in a direction opposite to rotor rotation. In an 
absolute frame of reference this speed of circumferential propagation 
is 30-60 per cent of the rotor speed [Ref. 8]. When these conditions 
become severe enough to cause net flow-rate fluctuations with time, the 
Situation is called "surge." This describes operation in the unsteady 
area of the compressor performance map in Fig. 1. Valensi [Ref. 9] has 
made a thorough investigation into these phenomena with some excellent 
photographic results. Of the various analytical methods available for 
predicting compressor stall, the most generally accepted employs the 
limiting-blade-loading criteria of Lieblein, Schwenk, and Broderick 
[Ref. 10]. Their "limiting diffusion factor," D, is a function of the 
velocity relative to the blade and the circulation about the element. 
If a constant axial velocity assumption is made, this simplifies to a 
function involving the flow turning angle and cascade solidity. Рог 
design purposes where such stage parameters may not be available, 
Lieblein [Ref. 11] has derived an equivalent diffusion ratio, Deg? which 
is based on inlet and outlet conditions to the compressor and the type 
of blading installed. Unlike a single airfoil where angle-of-attack is 
used to predict stall, the critical variable for a cascade of blades is 
the flow turning-angle. The exit-flow angle as well as the inlet angle 
must now be considered. 
2. Secondary Flows in Axial Compressors 

Returning now to порта compressor operation, flow patterns 
based on axisymmetric flow distribution and an infinite number of blades 
have been described quite well in the literature. However, it is necessary 


| 


| 
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to consider a finite number of blades since there are wall boundary- 
layer interactions, tip clearance leakages, blade boundary-layers, and 
centrifugal effects which are not represented adequately by assuming 
axisymmetric flow. Although these asymmetric phenomena have been 
studied independently for reasons of simplicity, it must be remembered 
that their effects are not necessarily additive (non-linear) and may 
have strong interactions. 

Figure 8 illustrates the various secondary flows which have 
been found to exist in axial-flow turbomachines as a result of the above 
mentioned viscous and three-dimensional effects. Lakshminarayana and 
Horlock [Ref. 12] have made athorough review of secondary flows in axial- 
flow turbomachines and have indicated where in the literature more 
specific information can be found. Woods [Ref. 13] presents the current 
state of the art in analytical treatment of secondary flows and associ- 
ated losses in axial-flow turbomachines. 

The induced velocities caused by secondary flows are considered 
losses since only a small part of their kinetic energy is recoverable 
downstream. Secondary velocities may perturb the exit flow angle from 
a blade row. Near annulus walls where lies the origin of the central 
vortex in Fig. 8, the induced exit flow angle can cause local blade 
stall. Much work has been performed in analyzing these phenomena, and 
some analytical methods for predicting three-dimensional exit flow 
angles have been proved valid through experimental results [Refs. 14 
and 15]. 

Another source of strong secondary flows in blade passages lies 
in non-uniform inlet flow such as described previously. Treating inlet 


flow distortion as vorticity provides a natural basis from which to 
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FIG. 8 | SECONDARY FLOW AND VORTICES IN AN AXIAL 
FLOW COMPRESSOR ROTOR [adapted from Ref. И | 


compute resultant secondary circulation at the exit plane. The effects 
of this secondary flow on cascade performance can be studied analogously 
to boundary-layer induced circulation. Comparison with an appropriate 
blade-loading criterion may predict the compressor's unstable operation. 
3. The Theory of Vorticity Analysis 

The thermodynamic properties which are significant in compressor 
performance are total-pressure (velocity), total-temperature (enthalpy), 
specific heat, and molecular weight. Variation of any one or combination 


of these properties with time and/or position at the compressor face 
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constitutes inlet flow distortion and should be an integral part of any 
universal distortion parameter. To date distortion factors dealing only 
with one inlet property at a time have been reported, i.e. total- 
temperature variation ^ Ky > Kr(stall) or total-pressure fluctuation > 
Kp > Kp(stall). This type of piecemeal approach, while providing satis- 
factory results for certain restricted applications, cannot be expected 
to remain valid under all conditions. In reality, an aircraft engine 
will encounter all of these fluid distortion properties at one time or 
another. A partial list of causes of these is: inlet characteristics, 
environmental conditions, rocket exhaust from ordnance, hot runway opera- 
tion, steam catapult launches, and the normal operating envelope of the 
aircraft. There is definitely need for a more realistic or universal 
distortion parameter, one which accounts for the total inlet distortion. 
It is the intent of this work to show that a distortion parameter based 
on vorticity would be more flexible and yield more comprehensive surge 
prediction criteria than present methods. 

Crocco's Theorem provides a versatile means of converting exist- 
ing inlet distortion maps to vorticity in that pressure, temperature, 
and time are all available as desired for inclusion in the transformation. 
As indicated previously, time is an important parameter in such an analy- 
515, the aspects of steady-state versus turbulence induced stalls being 
quite distinct. Hawthorne [Refs. 16 and 17], Preston [Ref. 18], and 
Horlock [Ref. 15] all deal with the passage of vorticity through a cas- 
cade. Horlock and Hawthorne calculate induced exit flow angles which 
involves a solution to Poisson's Equation. Calculating appropriate exit 
circulation and comparing it with blade-lift circulation is a far simpler 


process that yields a qualitative feel for the relative size of this 
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effect. Comparison of this secondary circulation with that caused by 
three-dimensional effects should also be made in order to look for 
possible interactions. 

The effect of inlet distortion on blade angle-of-attack must 
also be considered since this is the other "half" of the flow turning 
angle. This effect is much less subtle and can be dealt with in a more 
straight-forward manner. Areas of low pressure represent lower axial 
velocity, hence higher angles-of-attack. The unsteady response of the 
/ 


blades must be considered, particularly for circumferential inlet 


distortion where blade passing frequency is a factor. 
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III. CONVERSION OF PRESSURE TO VORTICITY MAPS 


A. CROCCO'S THEOREM 


In general form, Crocco's Theorem can be written as 
T grad S + u x curl u = grad H + = (1) 


which holds in unsteady flow [Ref. 19]. From basic Thermodynamics 


5-5 р Т 
EN Со. 2 ==) E (=) (2) 
Y - 1 


ў "n р m 


Combine Eqns. (1) and (2), assuming H = Constant 


-RT 3u 
EB grad P+ t u x curl u = = 3 
Pr 9 quee 3 (3) 


The following development will apply to the system defined in Fig. 9. 


0 
202 3 Ы IDEE. ы 
peo oce 

im Vtip 
Shock —/ 
Wave — K—— — L —— М- No. of Blades 

С = Свога 

FIG. 9 |. NOMENCLATURE 
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Writing Eqn. (3) in component form using cylindrical coordinates yields 


-RT 9PT ~ 2 Чу a 

E ED e, + (чеш, = ЛЭГ = p er (4) 
T 

-RT 1 Py A A ди. 2 

En Lower ur ” 

-КТ ЭРЧ ^ ^ dU, A 

тата 2 о 7 н 5 2 


The following reference quantities taken at the compressor face will 


be used to non-dimensionalize the equations. Refer to Fig. 9. 


Р = Б. Pr = Average Stagnation Pressure 
u= uta a = Average Speed of Sound 
r= “2 D = Compressor Diameter 
ү 
Lo 2 
ш =! E w = Fluid Vorticity 


Table I lists some of the characteristic times available. Blade passing 
frequency will be used in this development, since most available data 


include compressor speed; then 


ў 2-5 21 


МИРУ 


В 


The numerical value of the reduced frequency for stall is based on the 


work of Carta [Ref. 20]. For values less than approximately 0.1, blade 
| 


29 





TABLE І CHARACTERISTIC TIMES 


BLADE PASSING FREQUENCY 


VripN (500) (36) 











= ————— = 4500 Н 
2 т (6) (2). р 
ONCE-PER - REV 
М 1500 
MP - = 
228 (62) ~ 125 НЕ 
WAVE PROPAGATION 
a І200 
RADIAL ма. ЭГ. 
OR 4 ‚300 Hz 
atu 
LONGITUDINAL [- * 60,80 = 120 Hz 
STALL TIME 
REDUCED | тіс 20, 
FREQUENCY U | 


д (0.1) (1200) _ 
| (3) (0.2) — 
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stall is essentially equivalent to steady-state stall. For a reduced 
frequency larger than 0.1, stall is suppressed due to dynamic effects. 
Other work in this area has been performed by Carta [Ref. 22]. The 
paper by LE and Reding [Ref. 21] is a survey of dynamic stall. 
Rewriting the component equations in non-dimensional form and 


canceling terms yields 











Эр" 22 Ve. Na узи ' 
1 а ТІР | г.) 
-— — + —  (u,'w_' - ujw,') = —— | — (7) 
ХР' дг' а2 9 2 = 2та2 ә%! 
1 а pud ! 
1 ЭР а УттрМа ди 
- БИН e cj e (25 (8) 
"Du 99 42 га Zo 2та2 ot! 
! = Pe ! 
1 ЭР а V Ма /ди 
1 1 1 1 2 
-----%-; (шыт - аушы) = | | (9) 
ҮР" 32" а 2та 2%! 


Equations (7), (8), and (9) are exact. Orders of magnitude of various 
terms will be examined to ascertain whether or not complexity can be 


reduced in a logical and consistent manner. 


B. ORDER OF MAGNITUDE ANALYSIS 
Figure Д-2© represents an inlet-flow distortion map at the face of 
a J85-GE-13 turbojet. This distortion map triggered surge. Тһе пар 


of Fig. A-3 is a similar situation in a TF-30 turbofan engine. These 


Reference to figures and equations labeled with a letter/number 
combination have more detailed explanations which are contained in the 
Appendix. 
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maps will be considered representative of stall-inducing flow distortion 
for purposes of making an order-of-magnitude analysis on Eqns. (7), (8), 
and (9). 


Establish an order of magnitude, € , such that 


шэг) 


й! 
From Figs. A-2 and A-3 
020 


Between 12 and 5 o'clock, Fig. A-2, 3P'/ar' is a maximum. Refer to 
Appendix A for all numerical values. 


ӘР _ a(P-/P-) _ ШЕ 524997 


^ b 2 (2r/D (7.1 - 6.05)(2/16.1) 


I 
= 


293 


From Fig. A-3 at the 12:30 o'clock position 





әр! 1.04 - .96 
— 2 = .685 
ar” (2.15 - 1.8)(2/6) 
ӘР! 
-0.7 < ul. > 
ar 


Between 5 and 7 o'clock, Fig. A-2, a maximum value of 3P'/30 exists such 


that әр" 1.20 - .95 
А -- = = .416 per radian 
90 .6 rad 


Similarly from Fig. A-3 at 1 o'clock 


ӨР! 1.04 - .92 
— = —————— = .343 per radian 
99 .39 Kad 

EN a 

99 
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A discussion of the term 3P'/83z' is contained in Appendix A, Part II. 


An acoustic analysis taken at the compressor face yields Eqn. (A-2) 


which allows a magnitude estimation of əP'/əz'. 


From Figs. A-3 
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through A-8 at the point "A", the following data are 


= .92 Pr = P'P- = 1480 psf 
.96 PY = 1545 psf 
- 1/2 Е 1/2 
MOS) _ (1610 ~ 1440) 
р/2 .0015/2 
475 ft/sec 


и/М = 475/.4 = 1190 ft/sec 


3 ТҰ T | шин = - 1480 lbs/ft^ ` 
1610 lbs/ft? | 


ална p Ze C NE 
475 + 1190 6 Fig x .2x10 ~ sec/Fig 
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An average value of ЈЕ for the J85 can be obtained from data in 
Appendix С. 
и 460 


"ее 9.415 
а 1110 


и ' 
2 


A maximum value would be 


1/2 
P' = 1.17 uc (1:55) - 750 ft/sec 
о/2 
750 
Нас = S - .675 
US se 7 


Data with which to estimate the magnitudes of u,’ and u,' or their 
partial derivatives are not available. However, an upper bound of order 
e can be definitely established by considering the consequences of Un 


eeu. = 0(=0). A radial velocity component of order unity is not 


9 
possible within a duct. Any u,' components present at the compressor 
face are caused solely by turbulent fluctuations and are probably of 
order c^; however an upper bound от е will be assumed. A velocity 
component such as u,', called pre-rotation, is often generated intention- 
ally as a means of establishing a desired angle-of-attack on the first 
blade row. The magnitude of this velocity should be smaller than u,' 

and certainly not larger. A velocity in the circumferential direction 


0 


of order e”, about equal to the axial velocity, is unrealistic; and if 


1t were present, the engine would just not run. For the case under 


4 


consideration there is no pre-rotation so that any ug' which does exist 
is caused also by turbulent fluctuations and is probably of order ce. 


Again an upper bound of e will be assumed. 
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Therefore 


Partial derivatives of u¿' and М with respect to time can also be 


given only a bound with respect to и, 


i i i 
OU, ди. 7 


at' ot! Ор: 








There are two phenomena causing a time dependency. One is due to the 
motion of the rotor. For patterns of large extent on the Py Map, the 
appropriate frequency is once-per-rev. For patterns of limited extent, 
blade passing frequency is appropriate. Unsteady aspects of the flow 
cause a time dependency also. 


Continuing, the magnitude of ur 
wD ди д D 

EP E (2 24 

га 90 92 / га 


From Fig. A-2 between 8 and 9 o'clock, and data from Appendix C 


will be assessed. 


(? р yis 

S 
P' = ‚97 ш = \ - 400 ft/sec 
р/2 


Ри. 17 u 


И 


750 ft/sec 
48 = 30° = „52 radians 


(350 ft/sec)(1.34 ft) 
0. 5 --------------------- = .814 per radian 
(.52 rad)(.605 ft)2(1110 ft/sec) 
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Similarly, approximating w,' yields 


| | os D | au дц, ) р 
? 2а 


0 M 
9 2a ar 


From Fig. A-2 between 12 and 5 o'clock 


-(750 - 400) ft/sec 1.34 ft 
Qa = s+) = -3.52 


адбі е 2(1110 ft/sec) 


ШЕ - 25 


For the axial component of vorticity 


= - - | “е. -— = 


07 ETE = 
2a or 90 2а 


Since both partial derivatives are really unkown in size, this term 


will be evaluated later. 


The magnitude of ди, '/9Е' in Eqn. (9) can be estimated by again 


referring to the point "B" in Figs. A-3 through A-8. 


: 1/2 
(Pr - Р) V2 ((2)060 - 1440) 
P' - .92 ПЕЕ не = = | ----------- 
p/2 .0015/2 
uz = 231 ft/sec 
Р' = .96 US ОТЕ сес 


At = (6 figs.) x (.2х1073 sec/fig.) 
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ou," т E Ema 3/4 - 23] 
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Entropy can be expressed as 


5-5 H o. [T 
RENE тот 


where P and T are static values at the compressor face. Assuming P = P, 


combination of this equation with Eqn. (2) yields 


yl ца 
wc (PY = 1+ (6) 
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го 
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DE = lustai 


The deviation of local speed of sound from average speed of sound is less 
than the deviation of local stagnation pressure from average stagnation 
pressure. Therefore a/a will be set equal to 1. 

Table II is a summary of these established relative magnitudes in 
terms of е. The yelocities, Ug’ and u,', are placed at their absolute 


upper bounds and w,' is set equal to E for later verification. 


— —— ae + 


TABLE TI RELATIVE ORDERS OF MAGNITUDE 
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Now Eqn. (7) can be written 





1 ӘР! V N 3u,' 
— — | (: p - 1; ug!) ПРВИ не (7а) 
ҮР" ay! 9 21 а ot! 
ед ед € =0 =0 =0 =0 сё 
Dropping terms of order e 
1 әр! 
a. 7 - u,'u' = 0 
! ! 2- 9 
ҮР» Эг l (10) 
-] ЭР" 
ш, = БЕРЕ 
> yP'u 455) 
Similarly for Eqn. (8) 
1 әр! ү N 3u,' 
4 2 р (uz'u,.' 5 Up o, ' = Ца ш (8а) 
nn pe 27 а 9%! 
=0 е0 гд EU E =0 =0 c? 


1 ЭР! 
o A au 
yP'u,'r'\39 


The size of w_ was assumed to be of order 2 although at this point it 
is, in fact, unknown. In arriving at Eqns. (10) and (11), QM MS 


multiplied by small terms (u,' and ue n therefore little if any error 


r 
should be introduced here. Subsequent analysis establishes the size of 


w,'. From Eqn. (9) 





1 ЭР" ү М 54." 
Ee a a (12) 
үр! 92! Zn a gt! 
ы | | | | | | 
е0 € Е =0 € =0 0 € 
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Equation (12) contains most of the terms whose magnitudes were generally 
uncertain. All terms are of the order e so no consistent simplification 
can be made. This equation will not be used in the subsequent trans- 


formations from Рт to vorticity. 


C. APPLICATION OF RESULTS 

Equations (10) and (11) are the required relations for transforming 
total-pressure maps to circumferential and radial] vorticity maps 
respectively. The procedure, as applied to the J85 map, is outlined 
below. 

The expressions 3P'/3r; in Eqn. (10) and 3P'/936 in Eqn. (11) аге 
obtained in a similar fashion from Fig. A-2. For one particular radial 
line, values of P' are plotted against r', the slope at any point being 
the local value of эР'/аг'. This is done for sufficient radials to 


adequately coyer the entire map. Local values of P' can be taken direct- 


; can be calculated (assuming 


ly from the Py map, and from these u 
constant static pressure). Similarly in Eqn. (11), for one discrete 
radius, values of P' at a position, 9, are taken at intervals through 
360 degrees. These values are plotted, the local slopes being numerical- 
ly equal to 3P'/30. Other values necessary to compute e,' can be taken 
directly from the map or readily computed. See Appendix B for the 
calculations outlined above. 

The results of these calculations, depicted in Figs. 10 and 11 as 


we and w,' contour maps, present a two-dimensional picture of conditions 
at the compressor face. Since V*(Vxw) = 0 , it follows from Gauss' 
Divergence Theorem that 


> + 


#/ (о • г) дА = /// (ум · о) ду = 0 
А y 
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Application of this theorem to the maps of wy and we' contours yields 
a quantitative estimate of the axial vorticity, шу, which is reflected 
in Table II. These calculations are contained also in Appendix B, and 
the resulting map is given as Fig. 12. 

Radial and circumferential vorticity maps for the TF-30 total- 
pressure data in Fig. A-10 have also been generated. They are Figs. 


13 and 14. 
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FIG. IO RADIAL VORTICITY CONTOURS, J85 
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FIG. II CIRCUMFERENTIAL VORTICITY CONTOURS, J85 
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FIG.13 RADIAL VORTICITY CONTOURS, ТЕ-30 


44 





EN 





FIG. 14 CIRCUMFERENTIAL VORTICITY 
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IV. ANALYSIS OF VORTICITY MAPS 


A. RADIAL VORTICITY INLET EFFECTS 

In vector notation radial vorticity is aligned parallel to local 
blade spans of the compressor, tw, pointing outward radially according 
to the right-hand rule. The vorticity T Fig. 10 is concentrated in 
two areas, -wp between 5 and 7 o'clock and te, at 8 o'clock. In 
analyzing the effects of these local concentrations it must be noted 
that four distinct situations can be represented by a vorticity map such 
as this. These are listed in Table III. Vortex-sheets are aligned par- 


allel with the flow and vortex filaments are convected with the flow. 





TABLE III Possible Combinations for Radial Inlet Vorticity 


Case Vortex . Blades Flow 
I Sheet Stator Steady 
I] Sheet Rotor Uns teady 
III Filament Stator Unsteady 
IV Filament Rotor Unsteady 
Case I 


A continuous sheet of vortex filaments passing through a stationary 
blade row (stator) is a steady-state problem. According to the Biot- 
Savart Law, velocities are induced in the fluid which cause changes in 


local blade incidences, Fig. 15. In à cascade r is the circulation 


about a blade. | 
| 
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Г = (М 1 - Vio) (13) 
and 
і = рү l (14) 
Тһеп 
AL ^p AV AT 
— = — + — + — (15) 
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FIG. 15 RADIAL VORTEX -SHEET PASSING 
THROUGH А STATOR 


Therefore blades in region "A" where AT/T is positive are more heavily 
loaded. 

The first stator row is behind a rotor so Figs. 10 and 13 cannot be 
applied here without knowing the influence of the rotor on w,.'. This 


r 
case could be applied directly to inlet guide-vanes however. 


Case II 


A steady vortex-sheet passing through a rotating blade row is an 


unsteady problem due to the relative motion between particle paths and 
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the rotor. Again, induced velocities cause variations in the absolute 
velocity, V1, which result in 481 (relative inlet flow angle) and hence 
AT and AL. MW is relative velocity in Fig. 16. As blade "A" passes 
through the vortex-sheet to position "B", it experiences a decrease in 
lift which is represented by a smaller Tıjjft. This results in shed 


circulation equal in direction to T and a resultant decrease in 


lift 
Strength of the trailing vortex-sheet. Dynamic response of the blading 
is a factor which must be considered in this case. 

For a rotating cascade, circulation associated with blade lift is a 


function of the peripheral velocity, M For the circuit abcd in Fig. 


16 this is 
Daft 7 5091 - Куә) 22 


The corresponding lift on this blade is 








С = ЛЕВАК = 05( 1 - M2 MW (17) 
Г үет 
W, 
B Wu 
i W 
Vi cre E 
V. 
v 4———— Shed Vortices 
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FIG? [6 RADIAL VORTEX-SHEET PASSING 
THROUGH А ROTOR 
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For our sign convention, -w 


causes blade unloading while +w, increases 


r r 


the local incidence. Comparison of circulation causing blade lift to 

that associated with inlet distortion (rgicg¢) from Fig. 10 should give 
a feel for the relative strength of this effect. For brevity, all mean 
quantities will be used and the dynamic response of the blading will be 


neglected. From data in Appendices C and D, for the J85 


Wa = 650 ft/sec 
Yu? = 321 ft/sec 
ыя 1.095 
Blade Gap (Sm) = a ft 
Blade Chord (C) = um ft 
АЕ = (И - М) (16) 
Е = 30 ft^/sec (Relative) 
gist ^ "rls (18) 
_ 2854.9 1 ЭР! 
o = (2) (0-1) 
P'/9.17P' -8.17 p 55.90 
From Fig. A-2 at 8 o'clock 
P avg A 
From Fig. 10 at 8 o'clock 
1 ЭР" 
тов --) = күз 
Г 90 
Therefore o, = + 372 sae (absolute) 
and Pane ne ft2/sec (absolute) 
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The ratio of circulations becomes 


pos 
шин - а 
dist 2.97 


which is about 10 per cent. It must be remembered that this is occurring 
at blade passing frequency, and any induced stall would be dynamic stall. 


Blades 


14,200 Rev\/31 Blades\_ 
хо си )= 1.3 = 


Blade Passing Frequency - ( 60 Sec Rev 


Case III 

Unsteady vortex filament motion through a stator is a more complex 
problem since time dependency of the fluid motion as well as the unsteady 
response of the blading must be considered. This situation can be 
pictured as the passage of a vortex filament or a local concentration of 
vorticity passing through a cascade. In a similar problem, Crimi and 
Reeves [Ref. 25] calculated sizable lift fluctuations as a helicopter 
rotor blade passed over a spanwise oriented vortex. Ward [Ref. 26] 
observed dynamic stall of a rotor blade which he asserted was caused by 
a radial vortex passing beneath the blade. Figure 17 depicts this 


Situation in a cascade. 


ДУ ee 
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FIG. 17 RADIAL VORTEX FILAMENTS 
PASSING THROUGH A STATOR 
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Case IV 
It follows that a vortex filament approaching a rotating blade row 
is also an unsteady problem. Figure 18 shows basically how W,7 is 


affected by wp. 





FIG. I8 RADIAL VORTEX FILAMENT 
APPROACHING А ROTOR 


B. CIRCUMFERENTIAL VORTICITY INLET EFFECTS 

Circumferential vorticity taken vectorially is positive in the 
clockwise direction. Due to the signs in Eqn. 10, the positive gradients 
in Figs. 11 and 14 represent -|wg'|. The circumferential vorticity map 
for the J85 in Fig. 1] shows a strong concentration of vorticity at the 
outer radii extending through about 300 degrees. Figure 14 shows a 
Similar characteristic for the TF-30 although it is weaker and less 
extensive. As in the case of w,', four basic situations can exist for 


possible analysis. These are listed in Table IV. 
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TABLE IV Possible Combinations for Circumferential Inlet Vorticity 


Case Vortex Blades Flow 
I Sheet Stator Steady 
II Sheet Rotor Steady 
III Filament Stator Unsteady 
IV Filament Rotor Unsteady 
Case І 


Circumferential vorticity is oriented perpendicular to the local 
blade spans of the compressor. A sheet of wy, approaching a stationary 


blade row is pictured in Fig. 19. 
“в 
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FIG. I9 CIRCUMFERENTIAL VORTEX- SHEET 
APPROACHING A STATOR 
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Induced velocities according to the Law of Biot-Savart cause V,] and 
hence oj to change ahead of the cascade. This results in variable lift 
along the blade spans which produces shed vorticity downstream. Shed 
vorticity will be dealt with in more detail later when considering 
secondary flows. 


Case II 


A circumferential vortex sheet approaching a rotating blade row 
presents two problems for analysis. Relative motion between the particle 
paths and the rotor gives rise to a streamwise component of vorticity as 


well as a component perpendicular to the blade spans, Fig. 20. 


dc 


W Rotor 
B Mug 
V Ч 
So 
FIG. 20 COMPONENTS OF CIRCUMFERENTIAL 


VORTEX- SHEET APPROACHING A ROTOR 


The concentration of wg at the outer radius in Fig. 11 is essentially 
constant for about 300 degrees of angular displacement. The life time 


of one particular characteristic such as this can be as long as the time 





for one revolution. The series of TF-30 РТ maps in Appendix A shows a 
life-time of close to 1/4 revolution. The significance of this is that 
one blade can experience the same velocity distortion for up to 360 
degrees of revolution which allows a steady-state analysis. Induced 
velocities caused by the normal component of this sheet are those shown 


in Fig. 19. Velocities induced by the streamwise component are pictured 


in Fig. 21. adii 
U Rotor 
и) 
5 
ENG. 21 STREAMWISE COMPONENT OF A 
CIRCUMFERENTIAL VORTEX- SHEET 
“АТ. А ROTOR 


The streamwise component causes the typical variations in lift along the 
blade span, but it also passes through the blade row and affects the 
exit flow angle. This will.be discussed later. The normal component 
also causes a spanwise variation in lift as shown in Fig. 19. 

For Wg of small circumferential extent, the Situation is unsteady 


due to relative motion between particle paths and the rotor. 
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Case III 

A circumferential vortex filament approaching a stator is an un- 
steady problem. The effect of the associated entry velocity perturbations 
shown in Fig. 22, a function of the dynamic response of the blading, would 
probably be negligible. However, as this filament passes through the 
blade passage, a streamwise component develops which can interact with 
existing secondary flows. This phenomenon will be discussed in the next 


section. 


TIP 


11) 


Ет» Induced Velocities 


РО се CIRCUMFERENTIAL VORTEX FILAMENT 
APPROACHING A | STATOR 


Case IV 

A circumferential vortex filament approaching a rotor is different 
than the previous case only because a streamwise component of the filament 
is present ahead of the blades. Again unsteady response of the blades is 


a factor to consider. 


ЭЭ 





C. PASSAGE AND EXIT EFFECTS 

Several authors have treated the problem of turning a non-uniform 
flow in a cascade of airfoils. Squire and Winter's method for calcu- 
lating secondary flows downstream of a cascade [Ref. 27] has been used 
extensively in the literature. Preston [Ref. 18], using a more simpli- 
fied approach, also treates the generation of streamwise vorticity which 
gives rise to secondary flow after a cascade. However, Preston assumes 
that in addition to the components of vorticity present in the original 
vortex filament, a change in circulation about the blade will cause 
circulation to be shed into the stream. Hawthorne [Ref. 16 and 17] 
attempts an analytical proof of the existence of these different com- 
ponents of vorticity in nonuniform flow, and he obtains reasonable 
agreement with experimental results. 

Consider an isolated vortex filament in a sheet normal to the flow 
direction as it passes through a cascade of stationary авои ее 
Fig. 23. Due to the difference in flow velocity over the suction and 
pressure surfaces of the blades, this filament will be turned and 
stretched as it passes through the cascade. Downstream of the cascade, 
segment ab' gives rise to a "distributed streamwise vorticity" through- 
out the passage exit plane. This description applies not only to a 
filament as shown here but also to a vortex-sheet. Squire and Winter's 


well known expression for this distributed vorticity is 


а aU 
ae га (19) 


where z is streamwise vorticity, e is the deflection angle, and Z is 


the spanwise position. 
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VORTEX FILAMENT PASSAGE 
THROUGH А CASCADE 


Segment b-b', part of the original filament a-b in Fig. 23, lies along 


the stagnation streamline and in considering the spanwise direction, this 


becomes a "vortex-sheet." This is called "trailing filament" vorticity. 


A third component of exit vorticity, "shed circulation" caused by vari- 


able spanwise lift due to the asymmetric inlet flow, also lies along the 


Stagnation streamline and gives strength to the trailing vortex-sheet. 


This component is labeled c-d in Fig. 23. 


Preston's simple analytical expressions for these secondary components 


of circulation, as applied to Fig. 23, are as follows. 


Distributed Secondary Circulation: 


dr, = E25 cos (7) dz 


And replacing £5 with Squire and Winter's expression for distributed 


secondary vorticity, Eqn. (19), yields 


dr 


1 


EX m cos (2) dz (20) 





dz 
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Trailing Filament Circulation: 


du 
Eo СЕ) dz 21 
2 (= ) ШЫ 2 ) 
Trailing Shed Circulation: 
du 
АЕ well 2-1 1 а» 22 
3 (в " ( 2 | (22) 


Trailing shed circulation is equal in magnitude and direction to dr. 
The net circulation associated with the trailing vorticity in one blade 
passage and from one airfoil in the form of a "vortex sheet" is the sum 


e 


of these three components. 


dr = у + dro + аг, 
аг = Ез [ 4 S sin (=) 2 оо 10 211 (23 
42 2 3 2 


From this exit secondary circulation the change in mean exit flow angle 
can be calculated. This total circulation is not uniformly distributed 
over the exit area [Ref. 18]. 

Investigations into changes in deflection due to cascade secondary 
flows have generally taken two lines of thought. Values of induced 
deflection at the lifting line were derived from the induced velocities 
calculated solely from the strength of the vortex sheets by Carter 
[Ref. 28] and Hausmann [Ref. 29] to mention a few. These methods neglect 
the effects due to distributed passage vorticity and hence are unlikely 
to provide satisfactory agreement with experiment [Ref. 12]. 

Hawthorne [Ref. 30], Smith [Ref. 31], and others have cons idered 
secondary motions to be induced by the distributed vortices within the 


passages. It is assumed that the flow downstream behaves as if each 


| 
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blade passage were separated from its neighbor by the blade wakes, i.e. 
vortex sheets. The secondary velocities in the blade passages, v 
circumferentially and w spanwise, are defined by a stream function, v, 


such that 


av - ay 
ду И 92 


The stream function satisfied continuity for this steady secondary motion 


superposed on the main flow. 


3y 92 
Substituting the stream function into the vorticity equation (uo = curl и) 
and setting it equal to Squire and Winter's expression for the distributed 


vorticity, Eqn. 19, yields Poisson's Equation. 


2 2 

9р р 95р _ ы du 
- - -&j - -2е 051 (24) 
әу ә22 2 а2 


Solutions to Eqn. (24) are tedious but have been accomplished by a series 
method [Refs. 27 and 30], by relaxation techniques [Ref. 17], and by 
electrical analogs [Ref. 32]. The latter reference considers variable 


spanwise deflection which causes Eqn. (24) to become 
Уф = -дє dz t Uy pre (25) 


An exit flow solution for a parabolic inlet profile similar to the case 


at hand is given by Honda [Ref. 33]. 
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Once secondary velocities, v and м, һауе been determined, local 
changes in the exit flow angle, Aap, can be computed. Vo is the main- 


stream exit velocity for uniform inlet flow with exit angle ao. 


hay = таг” (15) 


The exit flow angle for distorted inlet flow then becomes 


а = a5 + Aa 
dist Ё 


Comparison of 277 with appropriate blade loading criteria [Ref. 10 
and 11] could explain the compressor's unstable operation. 

Another method of dealing with secondary velocities, v and w, ob- 
tained from solving Poisson's Equation was used by Hawthorne [Ref. 17]. 
Again the assumption of blade passages extended downstream by vortex- 
sheets trailing the blades is made. There are no normal velocities at 
the vortex-sheets. The strength of these sheets can be determined by 
calculating Aw, the change in spanwise secondary velocity across them. 
As mentioned before, these sheets of fluid leaving the trailing edges 
are vortex-sheets composed of trailing shed and trailing filament 
vorticity; see Fig. 24. If the trailing filament vorticity is estimated 
and subtracted from the total strength of the sheet, the trailing shed , 


circulation remains. From this can be calculated the variation in lift 


on the blade since shed circulation and blade lift are related by 


ІШ С 
L pul’. 


Hawthorne's expression for this trailing filament circulation 15 


duj\f 1 sin 2a, - sin 2a 
dr, = (2) с мл 2 02001201. - 2є cos(a,) 26 
t dz E COS aj й (26) 


60 





cT Trailing Edge 
Y 





4--- 
| UROTOR 


FIG. 24 SECONDARY FLOW DOWNSTREAM OF A 
BLADE ROW [ adapted from Ref. 17] 


which is the same as Preston's formula given as Eqn. (21) except for 
different inlet/exit flow angles. Hawthorne's predicted shed circulation 
agreed quite well with an experimentally determined spanwise lift 


variation. 


D. CALCULATION OF SECONDARY CIRCULATION 
A calculation of these components of secondary circulation as caused 

by the circumferential vorticity in Fig. 11 was performed in order to 
quantitatively evaluate the effect of such distortion on a compressor. 
The analysis was made assuming steady-state conditions which means Case 
Il of Table IV, a circumferential vortex-sheet passing through a rotor. 
The following terms require definition. 

V - Absolute Velocity 

W - Relative Velocity 


0 - Rotor Peripheral Velocity 
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o - Absolute Fluid Vorticity 
с - Relative Fluid Vorticity 
2 - Rotor Angular Speed 
a - Absolute Flow Angle 


B - Relative Flow Angle 
> - Velocity Vector 
> - Vorticity Vector 


Subscripts; see Fig. 25 


r,@,z.- Absolute Coordinate Direction 
s,n - Natural or Streamline Coordinates 
S - Streamwise 
n - Normal to Streamlines 
0 - Upstream of Compressor 
1 - Compressor face (entry to blade row) 
2 - Blade-row exit plane 





Streamline 


FIG. 25 NOMENCLATURE 
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Horlock's method [Ref. 15] of treating the passage of a vortex 
filament through a rotating blade-row results in the following vector 
transformation. Figure 26(a) is the velocity diagram ahead of the rotor 
which is convecting the vortex filament, wy (one of many in a sheet). 
The absolute velocity is axial so that Wo, = -U орок: Figure 26(c) is 
the inlet absolute vorticity combined with rotor rotation to give the 
inlet relative vorticity. Equations (27) and (28) are the corresponding 
analytical expressions for these components of relative vorticity in 


terms of Gg. 


со] = +2 2 cos B4 - wg Sin & (27) 


222. сіп В] - % С05 8) (28) 


Тһе approach of a vortex filament to a rotating cascade and the resultant 
secondary relative vorticity at the passage exit are pictured in Fig. 27. 
This figure is similar to Fig. 23 except it applies to a rotor while Fig. 


23 is for a cascade. 


Combining inlet streamwise vorticity which passes through the cascade 
unaffected with the exit vorticity of Eqn. (19) yields the total dis- 


tributed vorticity in the blade passage exit plane. 
662 = 65] -2e 23 (29) 
Substituting Eqns. (27) апа (28) 


сер = (20 со5 8] - o, sin 84) - 2e(-w, соѕ Ву - 20 sin 81) (30) 


9 


Since this exit vorticity is in the direction of the flow, multiplying 
it by the blade passage exit area normal to the direction of flow yields 


exit distributed circulation. 
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FIG. 26 VECTOR TRANSFORMATION FROM 
ABSOLUTE TO RELATIVE VORTICITY 
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BIG. 27 PASSAGE OF A VORTEX FILAMENT 
THROUGH A ROTATING CASCADE 


dry 


dr ` 


Эр S cos 8; (31) 


When transformed from absolute to relative coordinates, Hawthorne's 


expression for trailing filament circulation given previously as Eqn. 


(26) becomes 


аг; 


1 $1п 28, = sın 28 
-— = 5 Becher |. 22 2€ cos B 32 
dr q (> COS By - me 


Assuming two-dimensional flow, Hawthorne's equation for shed circulation 
is 

ага 

а= СОБР) Qj (tan 85 - tan 84) (35) 


A simplified application of these equations to the circumferential 


vorticity and total-pressure maps for the J85 (Figs. 11 and A-2) was 
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made by assuming an average radial distribution of 3P'/ar' and P' to 
exist for a circumferential segment of the face. Equation D-2 was 
applied to radials at the three circumferential positions shown in 


Fig. 28. (Appendix D) 


We = Sr sec”! (D-2) 
P 


Попав туо 


Now using Eqns. (27), (28), and (29), the various components of relative 
yorticity were calculated for a sufficient number of radial positions; 
see Appendix E. Figure 29 shows the relative sizes of these components 


for circumferential position I. 





FIG.28 DATA POINTS FOR CALCULATION OF М, 
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Secondary circulations as a function of span at the blade row exit plane 
were calculated using Eqns. (31), (32), and (33); see Appendix E. The 
area under curve "A" in Fig. 30 represents distributed secondary circu- 
lation. The area under curve "B" is trailing filament circulation, and 
the area under "C" is shed circulation. These latter two are the tir- 
culations due to the vortex-sheet. The magnitudes of these components 


are found by integration over dr'. 


.9 
R Í t2 аг = 32.1 ft2/sec 
3 dr 


where R is the radius used for non-dimensionalizing. 

Concerning the results in Table V, distributed secondary circulation, 
T1, represents à loss in kinetic energy. Blade-lift circulation is 
approximately equal to 30 ft2/sec; therefore the ratio of these two is 
about 1 to 1. In order to evaluate the actual effect г. has on the 
compressor, a solution to Poisson's Equation, as described previously, 
would be necessary. From this solution the exact strength of the 
trailing vortex-sheet could be determined. Subtracting trailing fila- 
ment vorticity (or Го) from the calculated strength of this vortex-sheet 
would yield the exact shed circulation. Equation (33) from which Га 
in Table V has been calculated was found by Hawthorne to yield values 
higher than those found experimentally. However, this analysis is 
essentially an order-of-magnitude check, and Eqn. (33) is adequate. 


Since shed circulation is related to the spanwise variation in blade 


lift by 
AL * oW, г, 
іһел AL shed (34) 
east 
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The circulation about one blade in the rotor under consideration for 


uniform inlet flow is approximately 
Pipe = S(M4 - Wo) = 30 ft%/sec 


The average value of positive shed circulation over the three circumfer- 


ential data points considered, Radials I, II, and III, is 


авіш 3-9 151 
3 


8 ft2/sec 


So the variation in lift from Eqn. (34) is 


=e = 25 


which represents unloading at the center span of the blade. For the 


negative shed circulation at the blade tip 


-6.0 - 8.5 - 10.9 _ AL _ -8.8 
3 a 30 


If the compressor blading is operating near stall it certainly seems 


plausible that a change of 25 to 30 per cent in lift could trigger stall. 


/0 
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ү. SUMMARY AND REVIEW 


Figures 31 and 32 are presented as a review and eine of the 
preceding theory for analyzing the effect of inlet flow distortion on 
an axial-flow turbomachine. Data from total-pressure maps, 3PT/3r and 
3PT/30, consitute the initial known quantities. All pertinent assump- 
tions and simplifying restrictions are included in the chart as they 
apply. The resulting determination of w, through a consistent order- 
of-magnitude analysis implies complete definition of the velocity field 
ahead of the compressor with the exception of the ambiguity in sign of 
wz. This is significant in that a three-dimensional flow field has been 
determined from a two-dimensional pressure map at the expense of intro- 
ducing an error of order «2. Continuing in Fig. 32, the analysis next 
relates portions of several existing theories on secondary flows in 
cascades and turbomachines to the maps of radial and circumferential 
vorticity. This yields an estimate of the strength of the secondary 
flows caused by inlet flow distortion. The large value for AL/L which 
is calculated suggests that the linear cascade theories used are Dos 
adequate for the accuracy required in deriving a stall criterion. Some 
necessary considerations for further refinement of the theory which were 
not considered in this analysis are (1) Des twist, (2) three-dimensional 
effects and interaction between normal secondary flows with distortion 
induced flows, and (3) changes in vorticity with passage through the 


blade row. 
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Absolute Vorticity, wy Map 


Relative Vorticity Components: zy]; Es] 
Horlock [Ref. 15] 








552 = (т) 
Squire & Winter [Ref. 27] 


Total Exit Streamwise Vorticity, 552 


Circulation Components, Г Poisson's Eqn. > Induced Velocity 


Hawthorne [Ref. 16] y = Normal to Span 





Preston [Ref. 18] W 


Spanwise 





A) Vortex Sheet » AL Change Exit Flow 
AL/L = Ol hog/! Hawthorne [Ref. 17] Angle, Ag» 


ІШІ 32. 






Limiting Blade Loading 
Criteria, Lieblein[Refs. 11 & 12] 


UNIVERSAL DISTORTION FACTOR 


Composite Theory for Analyzing Vorticity Maps 
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VI. CONCLUSIONS 


This analysis is the first attempt at relating total-pressure 
distortion maps to the internal aerodynamics of a turbomachine. Sepa- 
rating the inlet distortion into components of radial and circumferential 
vorticity introduced flexibility in treating the problem as well as 
maintained a comprehensible physical picture of the aerodynamical 
situation. For the data analyzed herein, radial vorticity was found 
to be less significant to degradation of compressor performance than was 
circumferential yorticity since the radial vorticity was coupled to the 
dynamic response of the blading. This would require an unsteady 
analysis. The relative significance of Фу апа We in this analysis does 
not necessarily apply in general. The effect of axial vorticity, 972 
must be included in further work. 

The size of the calculated shed circulation from a blade indicates 
distortion induced secondary flows are of sufficient size to seriously 
affect the compressor's operation. With further refinement to the 
present theory, secondary flows could be evaluated to a degree of 
accuracy capable of predicting stall through the mechanism of shed 
circulation or mE in exit flow angles. A solution to Poisson's 
Equation would provide the needed accuracy with which to use appropriate 
blade-loading criteria in defining inlet-distortion stall criteria. 
Although stall normally originates somewhere downstream of the first few 
stages, the distortion must be treated as it passes through each indi- 
vidual stage. The present theory deals only with the first blade row; 
howeyer, it is one row closer to the origin of the stall and has elimi- 


nated a certain degree of uncertainty. This theory may be extended to 
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succeeding stages in the compressor by the technique illustrated in 
РГО. 26. 

Finally, in working toward a universal distortion factor, the 
effects of other types of inlet distortion, i.e. total-temperature and 
molecular weight, must be included. These extensions to the present 
theory should provide reliable and comprehensive stall-inducing 


distortion criteria. 
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APPENDIX A 


ORDER OF MAGNITUDE ANALYSIS 


Part I Preparation of Total-Pressure Data 


The following data and calculations apply to the experiments per- 


formed on a J85-GE-13 turbojet as described in Ref. 2, a = 0° stall 


point. Some data were also obtained from Ref. 24. 


subscripts; 


0 - 
1 
2 


Т 


, = 1380 psf 
Tr, = 532°R 
М, = 2.5 


Compressor Inlet 


Figure A-1 is a pressure contour 


Free stream 
Compressor face 
Rotor exit, first stage 


Stagnation quantity 


(Static) 


flow area, А] = 185 in? 


map of Pr1/ Pro reproduced from 


Ref. 2. For our purposes, this must be converted to P' which is defined 


as Рт1/Рт] : 
Pressure Recovery of Inlet = .77 

Р. 
E 7 
Pro 

Er. Im oe 

20 ТО Pn 

ШҰ Т 
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Application of Eqn. (A-1) to Fig. A-1 yields the desired pressure 
contours, Fig. А-2. Figures A-3 through A-11 are instantaneous pressure 
contours (P') taken from experiments conducted on a TF-30 engine [Ref. 


23]. The radial lines at the centers represent compressor rotation. 


RPM * 10,000 


1! 


time between Figs. ( sec )( a гч 
103 Кеу 32 5ед 


3 


ДЕЛЕ х 10 
= .2 ms. 
- Elapsed time (A-3 > A-12) = 1.5 ms. 


Figure A-11 is the frame prior to engine surge. 


Ре 1440 psf 


Estimated Inlet Values 


М = .4 


Ш 1 
— _ 5 = ==> A 
Brees В (1-5 яг) 


Р; = 1610 psf 
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Part II Pressure Fluctuations in the Z Direction 

In the case of steady flow, APT is present only when dissipation 
occurs. Some causes of such dissipation in an inlet are shock waves, 
flow separation, obstructions, and boundary layers. These examples 


are illustrated in Fig. A-12. 
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The stagnation pressure gradient is present only in the vicinity of dis- 
sipation. Therefore in steady flow, əP'/əz' is zero at the compressor 


face if the dissipation is located upstream as shown in Fig. A-13. 


м 


Streom Tube 





E16, A-13 PRESSURE HISTORY FOR 
STEADY INLET FLOW 
For unsteady flow, pressure gradients are again caused by dissipation. 
If there were a standing shock wave in the inlet, unsteady flow could 
cause this to oscillate in the axial direction. This is called "shock 


wave twitter" and can be self excited in oscillation. This is depicted 


in Fig. A-14. 
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Shock Wave 


FIG. A-14 SHOCK- WAVE ‘TWITTER 
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The changing geometry seen by this oscillating shock-wave causes pres- 
sure pulses to be generated as shown in Fig. А-15. А/А% 15 the area 
ratio as a function of x. The second quadrant is local Mach number as 
a function of area ratio, A/A*. The third quadrant is the stagnation 
pressure ratio for a shock at M. Finally, the fourth quadrant yields 


the derivative of Pr with respect to x. 


Д/д 





Pro 


FIG. A-I5 TOTAL PRESSURE FLUCTUATIONS 
FROM OSCILLATING SHOCK-WAVES 


The resulting АР /Ах, a function of diffuser design, can be evaluated 
at the compressor face by assuming an acoustic analysis with no wave 
reflections at the compressor. Figure A-16 is a plot of the pressure 
pulses (waves) in the x-t plane. The slope of the characteristics is 
equal to the reciprocal wave speed; curvature is due to the variable 


area at the center-body. 
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FIG. A-16 PRESSURE WAVES IN THE x-—+t PLANE 
AS A RESULT OF SHOCK-WAVE TWITTER 


Consider the three points labeled 1, 2, and 3 in Fig. A-16. 


ӘР" | 250 
920 - 224-241 
since Pz" = Р»! 
ap’ РЗ’ -Р' Pg! = Ру 
at! | 23 - t | 23 - “ 


equating 


93 





22' зан _ (2, - 21) (2/0) _ E (25 - 2, 
ар (#5 - 1) (2а/0) а (4-4) 
1 
а. = xg ШІ 
à 
therefore 
ЭР" 1 ЭР" 
ВА О) Site 
and 
ӘР 1 9Р 
СТ m и (A-2) 
92 (a + u) Sr 


The usefulness of Ғап. (А-2) сап һе demonstrated by application to a 
simple inlet/engine configuration with an oscillating shock wave present; 
see Fig. A-17. For brevity the following assumptions are made 

(1.) a = Constant 


(2.) Wave impedance of compressor matches wave impedance 
of duct ` 


(3.) U = constant 


(4.) -No wave dispersion 


Se c ^ 


— 
Shock ———* = 
Wave — и * 
*N m- 
> 
>= 
= | 
Хо ЖЕ 


FIG. A-17 NOMENCLATURE 
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The shock wave, oscillating within IN, produces pressure waves which 


are represented as a sine wave function. 


X-Xg = Nsin (ot) 
9х 
S = wN cos (wt) 


within AN 


123 9х је 


ot ot 9X 


д 
w N cos (wt) ES 
9X x 


0 


Using a wave transit time equal to L/(atu) with a = wave speed relative 


to the gas, the unsteady pressure at the compressor face becomes 








әк - v N cos u(t - о) 20 
ЭГ ати’ эх. 
Хүн, Хү 
Combining this with Eqn. (A-2) yields 
SD o N cos w(t - Ll) әр 
ч = atu 3 
Тасе m uta 2 


Non-dimensionalized the preceding equation becomes 


Шоо: | Le 


92! u 9X 


where 2 is rotor angular speed. 
This is one method of relating pressure fluctuations at the compressor 


face (3P'/az') with the pressure waves generated by inlet/shock-wave 


interaction. 
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APPENDIX B 
CONVERSION OF Py MAP TO VORTICITY - 485 


Pressure data taken directly from Fig. B-1 are listed in Table B-I. 
From these data the curves in Fig. B-2 are generated. Slopes taken 
from these curves, which are tabulated in Table B-11, provide a means 
for constructing the circumferential vorticity map in Fig. 1l. A 
similar procedure is performed in generating the radial vorticity map 


in Fig. 10. Data are contained in Tables B-III and B-IV and Fig. B-3. 
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PRESSURE DATA AND SCALES FOR J85 


AiG, Bal 
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Table B-I Pressure Data Taken in Radial Direction from J85 


Total-Pressure Map 
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.92 1.04 
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Table B- II 


09 
op'/ar' 
0 

.09 
. 28 
2 
1.8 
2.6 
1.6 
‚83 
.37 
0 
-.6 
22.7 
-3.0 


Pressure Gradients Obtained from Slopes in Fig. B-2 


30? 

r' əP'/ər' 
.4 16 
-5 .16 
.6 ‚26 
E 2b 
.76 1.4 
522: 
58 6 
222911: 
286 0 
.9 -.72 
.94 -2.7 
.96 -3 


60° 

r' oP'/or' 
4 207 
D e 
.6 188 
27 125 
i 55 
4. 223 
202. 425 
(db — 55 
.89 0 
.92 -.6 
.94 -2.7 
.96 -3.0 
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90° 

t aP'/ar' 
22 МБ 
0 ІМ 
.44  .26 
32 0 

.6  -.26 
‚68 -.27 
.66 0 

“Я . 29 
JE 122 
Ae 2.0 
Эй 2.0 
-88 16.0 
88  .9 
цо xe 


120? 

ис n 
.4 .03 
12 .06 
16 al 
35515 233 
“128 1.88 
.8 3 
22225 
284 1.08 
286 .4 
289 0 
9  -.26 
92 -.6 
.94 -2.7 
.96 -3.0 


150° 

ro" Pp" 
гЗ; 25222 
22 7222 
5 -.22 
o B 
44 .3 
.6 253 
21 d 
25 15 
.48 4.0 
.8 7 
ады 3 
225 0 
9° 280772 
.94 -2.7 





Table B-II (Continued) 
180? 210° 240° 2709 300? 380" 
нэг г Perm nt oP'7S$pr'L r'sP'/ar'l r' aP'/ar'| r' 9P'/or' 
E -.28 1.3 407 7.3 „1з 392 35057 .4 ‚35 | .4 ‚15 
‚4 "OSEE 3 “ЮУ... .13 | .4 Sur 25 (20810 233 
.48 -.12 [1.0 | | ТЕЙ 121127 708 
451 0 Ed MEC 0 SD M sz O 2462 1.9 
.6 «Эй 92 0 T sl. yg «е 5282 2. 
2/ 27 .94 -.15 | .76 1.6 TO 2.3 8 1.05 
Eoo 1.7 | 2/0. ^ 56 242173 400: “б 
| | a .44 | .82 .65 | .84 .23 
zo 1.7 3859: 0 О 06 0 
.88 0 (897 --23/:17:9: --,7 2 ee 
.9 -2.5 .92 -.6 9213 .94 -2.7 
.96 -1.2 .94 -2.7 .94 -2.7 .96 -3. 
.96 -3.0 96-3. 
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Table B-III 


0 
(rad) 
4.26 
6.17 


229 


Р! 


191 


23! 


T 


Pressure Data Taken in Circumferential Direction from J85 


0 
(rad) 
4.19 
4.48 
5 23 
.174 


.46 


Р! 


23 
197 
en 
221 


Total-Pressure Map 
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3.14 
4.10 
4.41 
5.48 


294 
.84 
. 84 
221 
297 
27 
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(rad) 
1.69 
2.41 


3:98 
4.36 


dte 


261 


Р! 


17! 
. 84 
.84 
91 
24 
297 


r' = .69 
0 p 

(rad) 

0 .97 
о 
2.79 .91 
3.84 .91 
4.31 .97 
4.59 1.04 
5.13 1.04 
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Application of Gauss' Divergence Theorem is performed in the follow- 
ing manner. A grid such as is shown in Fig. B-4 is placed over a 


yorticity map. 


d 





FIG. B-4 METHOD OF DIVERGENCE 
APPLIC ATION 
The net Wy passing through region abcd is equal to 





[(-.13)(83) + (759621 - TIGEN FAN > | 


Net flux of w¿' is calculated in a similar manner using the same grid 


over an wg map. The sum of flux of w,' and w¿' is then equal to 


w,'(ab-cd) so that 


юу (ар са) 
0 = Е 
$ Area (abcd) 


This procedure is performed on each of the 36 segments so that a complete 


wz' map can be drawn. These data are contained in Table B-V. If the sum 


108 





of net flux of ом" ала We is negative, then the net flux of a must 


be positive. There is an ambiguity of algebraic sign which cannot be 


resolved without additional data or information. One evaluates o'' n. 
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APPENDIX C 
ESTIMATION OF SOME J85 PARAMETERS 


= 
! 


Subscripts: Stagnation Quantity 

0 - Freestream ahead of Compressor 
1 - Compressor Face 

2 - Rotor Exit, first stage 


m - Mean radius 


Data obtained from Ref. 24: 


Pry = 1549 psf 
Pi = 1380 psf 
Тт] = 532°К 


Data obtained from Ref. 2 and Ref. 34: 


Т 
6 = 5 т (Corrected Total-Temperature) 
нэ ot (Corrected Total-Pressure) 
2116 pst 
N* = 16,500 RPM (Rated Engine Speed) 
и _ : 
Wcorr Татты” 15 KG/sec = 33.5 lbs/sec (Corrected Mass Flow Rate) 


Percentage Corrected Engine Speed = Цам х 100 - 66.9 


D4 = 16.1 inches (Compressor Face Outside Diameter) 
DI hub 7 4-8 inches 
A, - 155 іп2 (Entrance Flow Area) 
N = 31 (Number of blades in first row) 
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ТІР 


HUB 
Blades 
| 2? i ke 
j 24 tt 
5.24" 
805" 
FO CI BLADE PASSAGE DIMENSIONS 
Calculated Parameters 
P. a 
zn = (1 + 1-1 м2) ү-! 
Р. Y 
1 
1549 3.5 


ши - = 2 
1380 E 1.125 x (1 + 2M ) 


М -41 
— |- 
| 
| 
=. 
+ 
р 
=! 
= 
N 
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= 


As 


49.02 ~ 514 


o | 
i 


g | 
| 


1110 ft/sec 


- Ma 
us = (.414)(1110 ft/sec) 
U, = 460 ft/sec 
. 
RT} 2 
7 ft-lbs 
(1715.6 STug?R ) (514°R) 


p = _.001565 slugs/ft? 


a check on these figures, another method may be used. 





1549 ) 
ца а «аеш. (ззяр 1Бейзес) ( 2116 
S ,/532/518 
W = 24.15 Tbs/sec 


2 РА] 


(24.15 lbs/sec) (144 in@/ft2) 
s]uas паса 
(.001565 $1495 ) (155 112) (32.2 


СЕ 
N 
11 


462 ft/sec 


113 


155 ) 


Se ait | 


Slug 





Estimation of 81 





Eoo coe VELOCITIES AHEAD : 
OF ROTOR 


W > Relative Velocity 
V > Absolute Velocity 
A E EE 


а Уа 


Dy = .8/5 ft 


Ang. Speed = (86-9)(16500) 


14,200 RPM 
М Тт/518.7 X 100 


W = ОТБ 14200 геу 
ul т геу ) 60 гас 
Wad = 650 ft/sec Mean Radius 
в = tan? ( 298 ) Ву = 55° Mean Radius 
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Estimation of 82 


Individual stage characteristics and loss data for the J85 are not 
available. The following data are from Ref. 34: 

Compressor of Free-Vortex Design 

Lieblein Diffusion Parameter, D « .5 

Solidity (c) = 1.14 


ји Stage Rotor 
Axial Velocity Ratio = .96 


122). 0 (Ref. 10) 


Wy 2YW, 





Wo = Va2/cos 2 .96 V. 1/ C05 Bo 


= 442/cos Bo 546 


W 2 = Үә tan Bo .96 Ж tan Bo 


- 442 tan Bo 


802 ft/sec 


221 442/с05 82 E (650 - 442 tan 82) 
802 (2)(1.14)(802) 


Bo = 36° minimum 
Wig = 320 ft/sec 


Since the loading limit is D « .5, 36? can be considered a minimum for 
Bo, and the maximum turning angle becomes 
АВ - В; > Bo = 55° - 36° 
= 199 


€ - .3232 radians 


ІМ 





APPENDIX D 
CALCULATION OF ABSOLUTE VORTICITY 


Radial Vorticity: 


i 1 1 ӘР' 
ш ume m LS OF 11 
STEM bx "2 ў 
Y = 1.4 
Б 
prs E 
РТ 
| U U U 
2 = E = эн = [ES 
à Uz/M U, 
1/2 
m (error y [py - 1380 Y 
(PT-Ps)(2/p) 1549 - 1380 
He p 
M 1692 8517 
= Ја Е р 1/2 
: Ш BE. 8.17) - Н(9.17:1-8.7) 
- Pr 169 PT 


Ва. = MA/9.17P' - 8.17 


u) 
m 
a 


2 


Substituting into Eqn. (11) 


c 2a 1 oP' , ft/sec 
r 14DP' M4/9.17P:-8.17 Л 30 ft 


116 





Recall: ЕЕ ПОМО гс ёс 
16:1 
D = Ян ft 
M = .414 
2854.9 
и = Е) сес” (0-1) 
: Р' f9.1/P" - 8.16 \r' a6 
Circumferential Vorticity 
E -2854.9 ВА се! (0-2) 
Р" ./9.17Р" - 8.16 Хэг 


Equation (D-2) is used to convert the contours in Fig. 11 %0 
absolute vorticity for subsequent calculation of secondary circulation. 


Tabulated results are contained in Table D-I. 
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Absolute Vorticity for Three Circumferential Positions 


Table D-I 
© /r (in) 
23 2.43 
925 2.63 
235 2.83 
ЕЗІ: 3.03 
E 9223 
1225 3.43 
Б 3.63 
475 3:83 
5 4.03 
1525 4.23 
o NEU de 
2579 4.63 
„6 4.83 
(025 5.03 
Kos 5.23 
07/5 5.43 
Y 5.63 
„725 5.83 
.75 6.04 
.775 6.24 
.8 6.44 
625 6.64 
.85 6.84 
„875 7.04 
„9 7.24 
.925 7.44 
.95 7.65 
.975 7.85 
1.0 8.05 


p! 


on Wy” Map 


Radial I 


489 


895 
279 


| 
1. 
| 
1 


“IC Po Po © 


ES 
0 
-.25 
= 
Ip 
-1.5 
-3.0 


-2554. 
1200. 
-1567. 
-1416. 
-1297. 
-1201. 
-1121. 
-1053. 
-1053. 
-994. 
-994. 
-942. 
-942. 
1096. 
-855. 
-818. 
-1682. 
- 2156. 
- 2402. 
- 2243. 
- 1522. 
ІШЕ 


0 


>32, 
111125 
1522. 
2480. 
5608. 
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ONAP PP OTM MD NS OO NS SI OT OT PO PO PO ГО PR OY CO Co Ov CO *O 


aP'/ar' We (зес”!) E 


105 
¿05 
205 
.05 
.05 
‚06 
1075 
209 


Radial II 


-520 


-624. 
-/80. 
-656. 
“27. 
-802. 
-/67. 
-826. 
-/58. 
-860. 
-895. 
-1030. 
-1011. 
-1716. 
-2840. 
-3509. 
E 
-2283. 
-741. 


0 


%741. 
1522. 
5608. 
7179. 


9P'/oar' Wy (sec!) 


C OY MN г -J OY OY — N = WH C9 C9 XO iO O — — 


Со «о > ОУ 





Table D-I (continued) 


Radial 111 


Сы Ут си әР'/әг' ш (сес-1) 





3 2.43 .9 -.2  42080.3 
.325 2.63 

.35 2.83 

75 3.03 

75025 

‚425 3.43 

.45 3.63 

‚475 3.83 

‚5 4.03 -.2  42080.3 
„525 4.23 0 0 
55 4.43 105 
‚575 4.63 .4 -4160.7 
6 4.83 .5 -5200.8 
.625 5.03 ! .5  -5200.8 
ІШІ 52 ‚90 „5 

„675 5.43 90 5 

.7 5.63 9] .5 -3649.9 
.725 5.83 .94 1.0 -4478.9 
‚75 6.04 ‚97 1.5  -5149.8 
775 ООЛ 110 2.0 -5681.4 
8 6.44 | 1.04 2.3 -5380.8 
‚825 6.64 | 1.17 1.7 -2588.0 
85 6.84 11.18 0 0 
.875 7.04 | 1.2 “8-2 112522 
х 721 (115 205 124 
ЖЕТТІ e 1125  42283.6 
6 7.65 |104 -30 +7018.4 
975 7.85 | 1.0 -3.2  49090.3 
1.0 8.05 
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2.0 
_ ка 


БУ 


9 
474 
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APPENDIX E 
- CALCULATION OF CIRCULATION FROM og MAP FOR J85 


ts] = + 29 соѕ Ву - ug Sin &i (28) 
= Rev E і В 
+02 236.6 сес) cos 55° - wg sin 55 
ИІ Ct 74 - .8191 о, (Е-1) 
1 5 - 29 sin 84 + we cos 8; (29) 
sie 2) (deans “EV шый 55° + (cos 55°) w, 
зес 
сп] = - 387.6 - .573 и, В 52 
2 - 7-51: 2261) (29) 
- (271.4 - .8191 шө) * 2(.3316 rad)(387.6 * .573 шө) 
tgo = 528.4 - .439 и, | (E-3) 
аг 
--- 5 552 S cos 82 (31) 
dr 
1.095 in 
S = —————— s .09125 ft 
Ваау г С 
В» - 36? 
dry 
dn. = .0738 662 ft/sec (ЕЗІ 
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dr? m (1 510 282 - SimwéP m 2 (32) 
n 





dr e COS fi 
В] = 559 
Е = .3316 
аг? 
E ^ .0480 En] ft/sec (E-5) 
dr 
S 2 
m - S cos 8] £,, (tan 82 - tan 8) (33) 
dr, | 
= .0336 Zn] (E-6) 
dr 


Tabulated results for Eqns. (E-1) through E-6) using previously 


determined data are contained in Tables E-I, E-II, and E-III. 
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Table E-I 


"325 
‚35 
"375 
45 
2 
Ti 
ES 


.85 


295 
0 


Relative Vorticity and Circulation Components at Radial I 


! 
LO 
No) 
I 
E QQ сл го Го “М Со “М! “Y Сл Сл ГО ГО ГО ГО -Ь ОУ Со с) ОУ СО О 


551 


(ѕес- 1) 


+2364. 
Д 
11595; 
+143]. 
11531. 
+255. 
+1189. 
ToS. 
ма. 
+1085. 
+1085. 
+1043 
+1043. 
+1005. 
+972. 
+942. 
+1649. 
+2037. 
+2239. 
+2109. 
+1506. 
+1182. 
+271: 
-17. 
-639. 
-975. 
-1760. 
-4322. 


1 


> .ь ол о› © © 01 OY O о со 4> <> ~ ~ О О с с 4 4> 4> со 


7 
5 
1 
8 


on] 


(әес-!) 


+1076. 
+627. 
+5103 
23. 
+356. 
+300. 
+254 
+2 
+20: 
+182. 
+182. 
+152. 
1152. 
+126. 
1192. 
TOM. 
+570. 
+847. 
+988. 
+897. 
+484. 
+249. 
-387. 
-589. 
-1025. 
-1259, 
-1808. 
-3601. 


зь У «О СО СУ ОУ СО “М ХО Со ХО NINN В СООСООООӘО О О О С О Yu 


122 


SS 


+1650 


+1305. 
+1216. 
+1150. 
+1098. 
+1065. 
+1020. 
4990. 
4990. 
4964. 
4964. 
4942. 
+942. 
+922. 
+904. 
+887. 
+1267. 
+1474. 
+1582. 
todo. 
+1196. 
+1016. 
+528. 
+373. 
+40. 
-139. 
-560. 
-1933. 


O Y О о Су Јо ~ ~ ГО «О «О 00 O O — — бо 00 Y Y DW — —= O1 М! 


dr4/dr 
(secl) I(ft/sec) |(ft/sec) |(ft/sec) 


ТТЕ. 
+96. 
489. 
+94. 
+61. 
AW. 
HO. 
ЗЭ. 
+73. 
ша. 
ІШ; 
469. 
+69. 
+68. 
+66. 
+65. 
198. 

+108. 

+116. 
tu 
+88. 
+75; 
+38. 
TA 

22. 
-10. 
-41. 
- 142. 


“Y YY YY ЮО л О ФС Y СУ СО СО «л л O O1 01 P9 P9 — — CO XO О СО ~ о ~ 


Se 
29. 
24. 
au. 
TG. 
14. 
12: 


-171. 


лг СО СО Со СОО -— Со СО СУ Сә о + О Со СО го Го ОУ ОУ С) С) оор ХО —! сл Со Го 


dro/dr 


220: 
tede 
p 
+14. 
Tapas 
+ 19% 
+8. 
+72 
t 
+6. 
+0, 
57 
d 
+4. 
112 
A 
+ 
2200 
+33 
+30. 
+16. 
шэў 
-13. 
-19. 
-34. 
-42. 


-61 


lil 


с с» (п еп Со С? З (л (а (а “У P™ Сл ГО 


аг:/4г 


CO CO Q1 — O 0 MN -— CO 


а 
T 


14 


шигээ 
NO RO 


сл 





Table E-II 


285 


‚45 


290 


5 


who 


.85 


495 
1.0 


Relative Vorticity and Circulation Components at Radial II 


u) 


(sec^! ) 


-520 
-520 
-520 
-520 
-520 


-624. 
-/80. 
-656. 
-/29. 
-802. 
-/67. 
-826. 
-/58. 
-860. 
-895. 
1030. 
ШТ. 
1716. 
2840. 
3509. 
3738. 
2203. 
-741. 


ЛЕДІ: 
+1522. 
+5608. 
+7179. 


со О FO CO OY Oy r9 М — OY = ~ WOW WWM «О — 
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(зес”! 


+697. 


t2. 
+910. 
+809. 
7259. 
+970. 
+899. 
+948. 
HENO. 


1115, 
+1005. 
ЗИ. 
+1100. 
+1677. 
+2598. 
+3145. 
TOA, 
%2141. 

+878. 

+27]. 

-336. 

-975. 
-4322. 
-5609. 


) 


9 


сл Со сл СО - Со Со сл ЈУ О О ПГ CO X9 О о СО O nro -» co C9 


^n] 
(sec” 


-89. 


-29. 
+59, 
-11. 
+30. 
172 
+52. 
+00. 
+47. 
+105. 
11251 
4202. 
+192. 
+595. 
+1240. 
+1623. 
11 55. 
4920. 
+3], 

- 387. 
-812. 
- 1259. 
- 3601. 
-4501. 


D 
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Оз о (л бу бо «О = = Y O со с До Со Со + O) — Co МО 
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552 


(вес!) (ft/sec) I(ft/sec) |( ft/sec) 


1456. 


%802. 
+870. 
+816. 
+848. 
+880. 
+865. 
1891. 
+861. 
+905. 
%921. 
+980. 
+972. 
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